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DIRECT AND INVERSE RESULTS ON ROW SEQUENCES 
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Abstract. We give necessary and sufficient conditions for the conver- 
gence with geometric rate of the common denominators of simultaneous 
rational interpolants with a bounded number of poles. The conditions 
are expressed in terms of intrinsic properties of the system of functions 
used to build the approximants. Exact rates of convergence for these 
denominators and the simultaneous rational approximants are provided. 
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1. Introduction 

Let f = (/i , . . . , fd) be a system of d formal or convergent Taylor expan- 
sions about the origin; that is, for each k = 1, . . . , d, we have 

►^ ' oo 

r^ (1) fk(z) = ^2(pn,kZ n , <P„,k^C. 

n=0 

■^4- \ Let D = (Di, . . . ,Dd) be a system of domains such that, for each k = 

1, . . . , d, fk is meromorphic in D^. We say that the point £ is a pole of f in 
D of order r if there exists an index k £ {1, . . . , d} such that £ £ D^ and it 
is a pole of fk of order t, and for j ^ k either £ is a pole of fj of order less 
than or equal to r or £ Dj. When D = (D, . . . , D) we say that £ is a pole 
of f in D. 

Let -Ro(f) be the largest disk in which all the expansions fk, k = 1, . . . , d 
correspond to analytic functions. If i?o(f) = 0, we take D m (f) = 0, m £ Z + ; 
otherwise, -R m (f) is the radius of the largest disk D m (f) centered at the origin 
to which all the analytic elements (/fc,-Do(/fc)) can be extended so that f 
has at most m poles counting multiplicities. The disk D m (f) constitutes for 
systems of functions the analogue of the m-th disk of meromorphy defined 
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by J. Hadamard in [5] for d = 1. Moreover, in that case both definitions 
coincide. 

By Q m (f) we denote the monic polynomial whose zeros are the poles of 
f in D m ({) counting multiplicities. The set of distinct zeros of Q m (f) is 
denoted by V m (t). 

Definition 1.1. Let f = (/i, . . . , fd) be a system of d formal Taylor expan- 
sions as in (Jl}. Fix a multi-index m = (mi, . . . , md) € 7h% \ {0} where 
denotes the zero vector in Zl. S'ei |m| = mi + • • • + m<j. Then, for each 
n > maxjmi, . . . , md}, there exist polynomials Q, Pk, k = 1, . . . , d, such that 

a.l) degi-fc < n — m^, k = 1, . . . , d, degQ < |m|, Q ^ 0, 
a.2) Q(z)f k (z) - P k (z) = A k z n+1 + ■■■ . 

The vector rational function R nim = (Pi/Q, . . . , Pd/Q) is called an (n, m) 
Hermite-Pade approximation of f . 

This vector rational approximation, in general, is not uniquely determined 
and in the sequel we assume that given (n, m) one particular solution is 
taken. For that solution we write 

\Z) -t*-n,m — \^n,m,l) ■ • • i -ttn,m,d) — \±n,m,li ■ ■ ■ j P n,m,d) / ^n,rm 

where Q n ,m has no common zero simultaneously with all the P nt m t k and 
is normalized to be monic unless otherwise stated. Sequences {R n ,m} for 
which |m| remains fixed when n varies are called row sequences, and when 
|m| = 0(n), n — >■ oo, diagonal sequences. 

The study of simultaneous Hermite-Pade approximations of systems of 
functions has a long tradition (see [6]) and they have been subject to renewed 
interest in the recent past (see, for instance, [3] and the references therein). 
Many papers deal with diagonal sequences and their applications in different 
fields (number theory, random matrices, brownian motions, Toda lattices, 
to name a few). At the same time, few papers study row sequences. In this 
second direction a significant contribution is due to Graves-Morris/Saff in 
[8] where they prove an analogue of the Montessus de Ballore theorem which 
plays a central role in the classical theory of Pade approximation. See also 
[2]-[in] for different approaches to the same type of results as well as [TT] 
and references therein for least-squares versions. 

Before going into details let us briefly describe the scalar case (d = 1) 
corresponding to classical Pade approximation which is well understood. 
When d = 1 we write f = /, m = m G N, and R njm = R n , m - Given a 
compact set K C C, || • \\k denotes the sup norm on K. We summarize what 
we need in the following statement. 

Gonchar's Theorem. Let f be a formal Taylor expansion about the origin 
and fix m € N. Then, the following two assertions are equivalent. 

a) Ro(f) > and f has exactly m poles in D m (f) counting multiplici- 
ties. 
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b) There is a polynomial Q m of degree m, Q m (0) 7^ 0, such that the 
sequence of denominators {Q n ,m}n>m of the Pade approximations 
of f satisfies 

HmSUp \\Q m - Q n ,m\\ 1/n = < 1, 
n— >oo 

where || • || denotes the coefficient norm in the space of polynomials. 
Moreover, if either a) or b) takes place then Q m = Q m {f), 

(3) 9 = rjj) ' 

and 

l/n _ \\z\\k 



(4) limsup||/- J R r , 



<-n,m\\K n ,,,, 

where K is any compact subset of D m (f) \V m (f)- 

From this result it follows that if £ is a pole of / in D m (f) of order r, 
then for each e > 0, there exists no such that for n > uq, Q n ,m has exactly 
r zeros in {z : \z — £| < e}. We say that each pole of / in D m (f) attracts as 
many zeros of Q n ,m as its order when n tends to infinity. 

So stated Gonchar's Theorem does not appear in the literature and needs 
some comments. Under assumptions a), in [7] Montessus de Ballore proved 
that 

lim Q nm = Qm(f), lim R nm = f, 

n— »oo n— voo 

with uniform convergence on compact subsets of D m {f)\P m {f) in the second 
limit. In essence, Montessus proved that a) implies b) with Q m = Q m (f), 
showed that 6 < max{|£| : £ € V m (f)}/R m (f), and proved (HJ) with equality 
replaced by <. These are the so called direct statements of the theorem. The 
inverse statements, b) implies a), 9 > max{|£| : £ G T y m(f)}/Rm{f), and the 
inequality > in §S§ are immediate consequences of [H Theorem 1] . The study 
of inverse problems of Pade approximation was suggested by A. A. Gonchar 
in [U Subsection 12] where he presented some interesting conjectures. Some 
of them were solved in [12] and [13]. See pQ for a brief account of Gonchar's 
most recent results and a list of his publications. 

In [8], Graves-Morris and Saff proved an analogue of the direct part of 
Gonchar's Theorem for simultaneous approximation with the aid of the con- 
cept of polewise independence of a system of functions (for the definition, 
see [8] ) . They also established upper bounds for the convergence rates corre- 
sponding to ([3j) and (jlj). The Graves- Morris/Saff Theorem was refined and 
complemented in [2j Theorem 4.4] by weakening the assumption of pole- 
wise independence, improving the upper bound given in [8] for the rate ([3j), 
and giving the exact one for ([4j). Until now, results of inverse type for row 
sequences of Hermite-Pade approximants are not available. 
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Our purpose is to obtain an analogue of Gonchar's Theorem for simulta- 
neous Hermite-Pade approximants, characterizing the exact rates of conver- 
gence of the Q n:m and R n ,m- 

The underlying idea in inverse-type results is that a polynomial which 
is the limit of the denominators of the approximants must have as zeros 
the poles of the function being approximated, provided that the rate of 
convergence is geometric. However, the actual situation in simultaneous 
approximation may be rather complicated as the following example shows. 
Take f = C/1,/2), where 

n=0 n=0 

and m = (1, 1). It is clear that the unit circle is a natural boundary of def- 
inition for both functions f\ and fa and thus z = 2 cannot be a pole of f in 
any system of domains. However, results contained in [2] show that the de- 
nominators Q n ,m of the simultaneous Hermite-Pade approximants converge 
with geometric rate to the polynomial [z — 1/2) {z — 2). 

This kind of examples leads us to introduce the following concept which 
is actually inspired by the definition of polewise independence in [8] . 

For each r > 0, set D r = {z G C : \z\ < r}, T r = {z G C : \z\ = r}, and 
D r = {z G C : \z\ < r}. 

Definition 1.2. Given f = (/1, . . . , fd) and m = (mi, . . . , md) G % d \ {0} 
we say that £ G C \ {0} is a system pole of order t of f with respect to m if 
for each s = 1, . . . , r there exists at least one polynomial combination of the 
form 

d 
(6) ^Pkfk, degp k <m k , k = l,...,d, 

k=i 

which is analytic on a neighborhood of D\£\ except for a pole at z = £ of 
exact order s and there is no polynomial combination of the form ([6]) with 
those properties for s greater than t. If some component m^ equals zero the 
corresponding polynomial p\~ is taken identically equal to zero. 

The great advantage of this definition with respect to that of polewise 
independence is that we have liberated it from establishing a priori a region 
where the property should be verified. This turns out to be crucial. 

We wish to underline that if some component mt equals zero, that compo- 
nent places no restriction on Definition 11.11 and does not report any benefit 
in finding system poles; therefore, without loss of generality we can restrict 
our attention to multi-indices m G N , and we will do so in the sequel, 
except in reference to the convergence of the approximants themselves. 

Notice that the definition of system pole strongly depends on the multi- 
index m and that a system f cannot have more than |m| system poles with 
respect to m counting their order. During the proof of Theorem 11.31 below. 
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carried out in Section [3l we give a procedure for finding in a finite number 
of steps all the system poles of f with respect to a multi-index m under 
appropriate conditions. 

It is easy to see that a system pole may not be a pole of f or viceversa. 
For example, let f be the system given by ([5]) and m = (1,1). The point 
z = 2, which lies beyond the natural boundary of definition of f\ and fi is 
not a pole; however it is a system pole of f since j\ — fi has a pole at z = 2. 

On the other hand, take f = (/i, f2) with 

f - l + l f l 



1 z-2' J z-3' 

and m = (1,1). Then the points z = 1 and 2 = 3 are poles and system poles 
of f but z = 2 is only a pole because there is no way of eliminating the pole 
at z = 1 through linear combinations of f% and J2 without eliminating the 
pole at z = 2. 

To each system pole £ of f with respect to m we associate several char- 
acteristic values. Let r be the order of £ as a system pole of f . For each 
s = 1, . . . , r denote by r^ )S (f , m) the largest of all the numbers R s {g) (the 
radius of the largest disk containing at most s poles of g), where g is a poly- 
nomial combination of type © that is analytic on a neighborhood of Dui 
except for a pole at z = £ of order s. Then 

- R C, 5 ( f > m )= min ^,fc(f,m), 

fc=l,...,s " 

i?f (f, m) = Rf T (f, m) = min rv s (f, m). 

S = 1,...,T 

Obviously, if cZ = 1 and (f , m) = (/, m), system poles and poles in D m (f) 
coincide. Also, R^(f , m) = R m (f) for each pole £ of / in D m (f). 

By Qi m i (f , m) we denote the monic polynomial whose zeros are the system 
poles of f with respect to m taking account of their order. The set of distinct 
zeros of Qi m i(f, m) is denoted by V\ m \(f,m). 

The following theorem constitutes our main result. 

Theorem 1.3. Let f be a system of formal Taylor expansions as in ([T]) and 
fix a multi-index m 6 N . Then, the following two assertions are equivalent. 

a) -Ro(f) > and f has exactly |m| system poles with respect to m 
counting multiplicities. 

b) The sequence of denominators {Qn,m} n >\m\ °f simultaneous Pade 
approximations of f is uniquely determined for all sufficiently large 
n and there exists a polynomial Q\ m \ of degree |m|, Qi m i(0) ^ 0, 
such that 

limsup||Q, m |-Q n , m || 1 / n = 0<l. 
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Moreover, if either a) or b) takes place then Q\ m \ = Qi m i(f,m) and 

(7) fl = max{ jf :^GP| m |(f,m) 

If d = 1, Rn,m and Q n ,m are uniquely determined. Therefore, Theorem 1 1.31 
implies Gonchar's Theorem except for @ whose analogue will be presented 
in Section 3.2 to avoid introducing new notation at this stage. 

The paper is structured as follows. In Section [2] we continue with the 
study of incomplete Pade approximants initiated in [2] proving results of 
inverse type. Section [3] is dedicated to the proof of Theorem 11.31 and the 
analogue of (H|). 

2. Incomplete Pade approximants 
Let 

oo 

(8) f(z) = ^2<f> n z n , ieC, 

n=0 

denote a formal or convergent Taylor expansion about the origin. 

Definition 2.1. Let f denote a formal Taylor expansion as in (|8|). Fix 
m > m* > 1. Let n > m. We say that the rational function r n ^ m is 
an incomplete Pade approximation of type (n,m,m*) corresponding to f if 
r n,m is the quotient of any two polynomials p and q that verify 

b.l) degp < n — m*, degq < m, q^O, 
b.2) q(z)f(z) - p(z) = Az n+l + ■■■ . 

Notice that given (n,m,m*),n >m> m*, any of the Pade approximants 
Rn,m* , ■ ■ ■ , Rn,m can be regarded an incomplete Pade approximation of type 
(n,m,m*) of /. From Definition 11.11 and ([2]) it follows that R n ,m,k,k = 
1, . . . , d, is an incomplete Pade approximation of type (n, |m|, m/%) with re- 
spect to fk- 

In the sequel, for each n > m > m* , we choose one candidate. After 
canceling out common factors between q and p, we write r n ^ m = p n ,m/Qn,m, 
where, additionally, q n ,m is normalized to be monic. Suppose that q and p 
have a common zero at z = of order A n . From b.l)-b.2) it follows that 

b.3) degp n>m < n - m* - A n , degg n , m < m - A„, q n>m ^ 0, 

b.4) q n , m (z)f(z) - Pn , m {z) = Az n+1 - X - +■■■ . 

where A is, in general, a different constant from the one in b.2). 

The first difficulty encountered in dealing with inverse-type results is to 
justify in terms of the data that the formal series corresponds to an analytic 
element which does not reduce to a polynomial. In our aid comes the next 
result, which provides such information in terms of whether the zeros of the 
polynomials q n ^ m remain away or not from and/or oo as n grows. Let 

' n,m — \^n,l> ' * * ' ^>n,m„ J> n ■>_ 771, Tfl n ^ 777, 
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denote the collection of zeros of q n ^ m repeated according to their multiplicity, 
where degg n)m = m n . Put 



and 



S = sup inf {|Cn,fe| : n > N,m n > 1, 1 < k < m n } 

N>m 



G = inf sup {|Cn fc| : n > N, m n > 1, 1 < k < m n } 

N>m ' 



Finally, set 

r n = min{n - to* - X n - degp n , m , m - X n - m n }, n > to. 
From b.3) we know that r n > 0, n > to. 

Theorem 2.2. Let f be a formal power series as in (|8j). Fix m > to* > 1. 

The following assertions hold. 

i) // |A ra — A n _i| < to* — l,n > no, and S > then Ro(f) > 0. 

ii) If \(m n + X n + T n ) - (m„_i + A„_i + T n _i)| < m* - l,n > n , and 
G < oo t/ien either f is a polynomial or Ro(f) < oo. //, additionally, 
there exists a sequence of indices A suc/i i/iai degg nim > l,n 6 A, 
i/ien Ro(f) < oo. 

Proof. From definition 

(9) (g„, m / - P n,m) (Z) = Az n+1 ~ X - + • • • , 

and <? n ,m(0) / 0. 

We may suppose that inf {|Cn,fc| : n > Tio,m n > 1, 1 < k < m n } > and 
|A n — A n _i| < to* — 1, n > no- Normalize g raim as follows. If m ra > 1 take 



Win / \ 

\\ ( 1 - -z — ) = a nfi + a n> iz H h a n ,m n ^ mn , a njl 

fc=l ^ «»•*/ 



I •)]"[( 1 - -r— J = a nj0 + a n ,iz H h fl„, m „^ m " , a n , 1 

k= 

Otherwise q n , m ( z ) = 1 = a n,o- 

Using the Vieta formulas connecting the coefficients of a polynomial and 
its zeros it follows that there exists C\ > 1 such that 

(10) sup{|a ni fc| : < k < m n , n > uq} < C\ < oo. 

The coefficient corresponding to z k , k £ {n — m* — X n + 1, . . . , n — X n } in the 
left-hand side of ([9]) equals 

(11) 4>k + a n ,i(f>k-i -\ h a n ^ mn (pk-m n = 0, 

since degp njm < n — to* — X n . 

If m n > 1, (HO]) and CLE} imply that 

|0fc| < Ci(|0fe_i| H h |0 fe -m„|)- 

Therefore, for each A; G {n — to* — X n + 1, . . . ,n — A ra } there exists £;' G 
{fe — 1, . . . , k — to} (m n < to) such that 

(12) \<t>k\<C x m\(t>k'\. 
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Should m n = 0, for the same values of k, we have fa = and ([121) is 
trivially verified. Substituting n by n — 1, we deduce that for each k £ 
{n — m* — A n _i, . . . , n — A n _i — 1} there exists k' G {k — 1, . . . , k — m} such 
that 

(13) \fa\ < Cwilfa'l 
As n > no, we have 

n — A n _i > n — \ n — m* + 1 

and 

n — X n -i — m* < n — A n — 1, 

because |A ra — A n _i| < m* — 1. Consequently, the range of values taken by 
k due to relations (fl~2|) and (JT3J) are either contiguous or overlapping for 
n > riQ. Since n — \ n tends to oo as n goes to oo, we conclude that for all 
n>riQ there exists n' & {n — 1, . . . , n — m} such that 

(14) \<f>n\ < C im \ct> n/ \. 
Let A be a sequence of indices such that 

lirri \fa\ x ' n = limsup \fa\ x l n = l/i? (/)- 

"•SA n— >oo 



Choose n € A. Due to (|14p there exist indices n\ > ri2 > • • • > n Tn , n rn < 
no, where r n < n — no, such that 

\fa\ <C x m\(j) ni \ <•••< (Cim) r "|(/> nrn |. 

Consequently, 

l/i? (/) = lim \(f> n \^ n < limsup(Cim) r "/ n < C x m. 

neA n->oo 

Therefore, Ro(f) > (Cim.) -1 > 0, which proves i). 

As for ii), assume that sup{|C n .fe| : n > no,rn n > 1, 1 < k < m n } < oo 
and \(m n + A„ + r n ) - (m n _i + A„_i + T n _i)| < m* - 1, n > n . Set 
t„(z) = (z - l) r ". Define g n , m = t n q n , m and p n _ m = t n p n ^ m . Normalize q n , m 
as follows. If m n + r n > 1 take 

m n +T n 

Qn,m{Z) = [J (2 — Cn,/c) = ^,0^™" r ™ + ' ' ' + &n,m n +r„-l^ + b n<mn+Tn , 
k=l 

where b n fi = 1. Should m n + r n = we set q n ^ m = 1 = 6 nj o- Using the Vieta 
formulas, it follows that there exists C<i > 1 such that 

(15) sup{\b ri: k\ ■ < k < m ni n > no} < C2 < 00. 

The coefficient corresponding to z k , k G {n — m* — A n + 1, . . . , n — A n }, in 
the left-hand side of ([9]) equals 

(16) (/> fc _ mn _ Tn + b n ^4> k _ mn _ Tn+ i H h 6n,m n +r„0fc = 0, 

since degp,, m < n - m* - A n . 
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Should m n + T n > 1, (fT5j) and (|T6l) imply that 

\4>k-m n -r n \ < C 2 (\(l)k-m n -T n +l\ + ■ ■ ■ + \<t>k\), 

or, what is the same, for each k G {n — m* — A n — m ra — r n + 1, . . . , n — X n — 
"in — Tn}, we have 

|0fc| < C 2 {\4>k+l\ H h |0fe+ m „+T„l)- 

Therefore, for each A; € {n — m* — A n — m n — r n + 1, . . . , n — A n — m n — r n } 
there exists fc' € {A; + 1, . . . , k + nt} (m n + r n < m) such that 

(17) |0,| > JM. 

In case that m n + r n = we have </>& = for the same values of k and ()17p 
is also true. 

Using the assumption that |A n + m n + r n — A n _i — m n _i — r n _i | < m* — 1, 
it is easy to check, similarly to the previous case, that the range of values 
taken by the parameter k for consecutive values of n are either contiguous or 
overlapping. Also, n—\ n —m n —T n tends to oo as n goes to oo. Consequently, 
from (1170 we have that for all n > no there exists n' G {n + 1, . . . , n + m} 
such that 

(18) |*,|>JM 

Using (TT8l) we can find an increasing sequence of multi-indices {n s } se ^ + , 
n s+ i G {n s + 1, . . . , n s + m} and ni G {no, . . . , no + m} such that 

U I > l^i I 
'^ s+l1 - (C 2 m)*" 

Should / be a polynomial there is nothing to prove. Otherwise, changing 
the value of no if necessary, without loss of generality we can assume that 
<j) ni / 0. Then, 

liminf \6 n ., \ 1/ns+1 > , > 



limsup(C 2 m)' s/ns + 1 C 2 m' 

s—>oo 

since 

S S 

limsup < limsup = 1. 

s— >oo 1-s+l s— >-oo ni + S 

It follows that 

Rnif) = tj- < n < C 2 m < oo, 

limsup l^l 1 /™ " liminf |&>, +1 | 1/n ' +1 ~ " 

n— >oo s— >oo 

as we needed to prove. 

Finally, if / is a polynomial, say of degree N, we would have that for all 
n > N + m, f = p n ,m/Qn,m arid q n , m = 1- Consequently, if there exists A 
such that degg n)m > 1, n G A, / cannot be a polynomial and, therefore, only 
Ro(f) < oo is possible. □ 
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Lemma 2.3. A sufficient condition to have [A n — A n _i < m* — 1 and 
\(m n + A„ + T n ) - (m n _i + A„-i + r n -i)| < m* - 1 is that 

min{m n + r n , m n _i + r n _i} > m — m* + 1. 

Proof. In fact, for k = n — 1 and fc = n, if m/t + r^ > m — m* + 1 then 
< Afc < rn,* — 1 because Afc + m^ + r^ < m and the first inequality readily 
follows. On the other hand, 

\(m n + A n + T n ) - (m n _i + A„_i + r„_i)| 

= \(m n + A n + r n - m + m* - 1) - (m n _i + A n _i + T n _i - m + m* - 1)| 

and < TTik + Afc + Tfc — m + m* — 1 < m* — 1 for A; = n — 1 and k = n. 
Therefore, the second inequality also holds. □ 



Applied to Pade approximation (rn* = rn) , Theorem 12.21 and Lemma [ 
imply that if deg Q n ,m > 1 an d its zeros remain uniformly bounded away 
from and oo, for sufficiently large n, then < Ro(f) < oo. This result has 
not been stated elsewhere. 

Let us see some consequences of Theorem 12.21 and Lemma 12.31 on the 
extendability of a formal power series and the location of some of its poles 
in terms of the behavior of the zeros of the approximants. First we bring 
your attention to some results from [2]. 

Let Bbea subset of the complex plane C By U(B) we denote the class 
of all coverings of B by at most a numerable set of disks. Set 



o-(B) = M\^\Ui\ : {Ui}€U(B) 



U=l ) 

where \Ui\ stands for the radius of the disk £/j. The quantity o~(B) is called 
the 1-dimensional Hausdorff content of the set B. 

Let {(p n }n€N be a sequence of functions defined on a domain D C C and tp 
another function defined on D. We say that {<^ n }ngN converges in cr-content 
to the function <p on compact subsets of D if for each compact subset K of 
D and for each e > 0, we have 

lim a{z G K : \ip n (z) — <p(z)\ > e} = 0. 

n— >oc 

We denote this writing a- linin^oo cp n = <p inside D. 

We define the number i?^(/) as the radius of the largest disk centered at 
the origin on compact subsets of which the sequence {r n<m } n > m converges to 
/ in cr-content. In [2] we gave a formula to produce this number and showed 
that it depends on the specific sequence of incomplete Pade approximants 
considered. Set Z?^(/) = {z G C : \z\ < R^if)}- 

Among other direct- type results, we proved that 

(19) R m *{f)<R* m {f)<R m (f), 
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that Rm(f) > implies Ro(f) > 0, and that each pole of the function / in 
Dm(f) attracts, with geometric rate, at least as many zeros of q n ,m as its 
order (see [2j Theorem 3.5]). Therefore, Theorem 1 2 . 2 1 and Lemma 12.31 imply 

Corollary 2.4. Let f be a formal power series as in ([8]). Fix m > m* > 1. 
Assume that there exists a polynomial q m of degree greater than or equal to 
m - m* + 1, q m (0) ^ 0, such that lim^oo q Uym = q m . Then < Ro(f) < oo 
and the zeros of q m contain all the poles, counting multiplicities, that f has 
inD* m {f). 

We need a relaxed version of Corollary 12.41 for the proof of Theorem 11.31 

Lemma 2.5. Let f be a formal power series as in ([8]) that is not a poly- 
nomial. Fix to > to* > 1. Let r n ^ m = p n ,m/q~n.m be a n incomplete Fade 
approximant of type (n,m,m*) corresponding to f, where p n ^ m and q n , m are 
obtained from Definition \2.1\ and common factors between them are allowed. 
Assume that there exists a polynomial q~ m of degree m, q m (0) ^ 0, such that 
lim n _ i . 00 q n ^ m = q m . Then < Ro(f) < oo and the zeros of q m contain all 
the poles, counting multiplicities, that f has in Dj^(/). 

Proof. Let us show that the assumptions of Lemma 12.31 are verified for the 
incomplete approximant r n , m . Let r n , m = p n ,m/qn,m, where the polynomials 
Pn t m and q ntm are relatively prime. Since q m (0) ^ 0, then g n ,m(0) 7^ 0, n > 
no- Thus, p n ,m and q n ^ m do not have a common zero at z = and A n = 
for all n > no- As before, set m n = degq njm and 

r n = min {n-m* - degp njm , m - m n } , n > n . 

Notice that r ra = m — m n , n > no, because the polynomials q n ^ m and p n ,m 
are obtained eliminating possible common factors between q n , m and p n ,m 
and by assumption 

min {n-m* - degp n , m , m - deg <?„,„} = 0, n>n . 

Therefore, we have 

Tn n + T n = m > m — m* + 1, n > no, 

and Lemma 12.31 is applicable. 

From Theorem 12.21 we obtain < Ro(f) < oo. Now, from the fact that 
each pole of / in D^if) attracts as many zeros of q n ^ m as its order it follows 
that the zeros of q m contain all the poles, counting multiplicities, that / has 
mD* m (f). D 

In case that there exists R > R m * (/) inside of which / is meromorphic 
then Dji contains at least to* + 1 poles of / since D m * (/) is the largest disk 
where / is meromorphic with at most to* poles. We can prove the following 
inverse- type result. 

Theorem 2.6. Fix to > to* > 1. Let f be a formal power series as 
in ([8]) that is not a rational function with at most to* — 1 poles. Let 
r n ,m = Pn,m/q~n,m be an incomplete Pade approximant of type (n,m,m*) 
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corresponding to f, where p n ^ m and q n ^ m are obtained from Definition \2.1\ 
and common factors between them are allowed. Suppose that there exists a 
polynomial q m , of degree m, q m (0) / 0, such that 

(20) limsup ||g n , TO - q m \\ 1/n = 9 < 1. 

n— >oo 

Then, either f has exactly m* poles in D m *(f), which are zeros of q m count- 
ing multiplicities, or Ro(q m f) > R m *(f)- 

Proof. From Lemma l2.5l we have Ro(f) > 0. So, / is analytic in a neighbor- 
hood of z = 0. We also know that Ro(q~mf) > Rm*(f) since the zeros of q m 
contain all the poles that / has in D m * (/). Assume that Ro(q m f) = R m * (/). 
Let us show that then / has exactly m* poles in D m *{f). To the contrary, 
suppose that / has in D m *(f) at most m* — 1 poles. Then there exists a 
polynomial q m * , with deg q m * < m* , such that 

Ro{q-m*f) = Rm*{f) = Ro{qm*qmf)- 



Let 



then 



i0)/0) = ;a n z 



oo 
n=0 



R m *(f) = Ro(Qm*q m f) = 1/ limsup a/KI 



The n-th Taylor coefficient of q m * [q n ,mf — Pn,m] is equal to zero. Therefore, 
the n-th Taylor coefficients of q m *Qmf and q m *Qmf - q m *q n ,mf + q m *Pn,m 
coincide. Take < r < R m *{f) and recall that T r = {z € C : \z\ = r}. 
Hence 

-L / [qm*qmj qm*qn,mj \ qrn*Pn,m\y^) 



2m r UJ n+1 

•' 1 r 

1 f [Qm- Qn,m]{u) q m * (w)/(w) 



2ni J T oo n+l 

Making use of (|20p it readily follows that 

1 r . — r 9 

— — — ■ = limsup Vl a n| < -• 
Hm*{j) n— Xx> T 

Letting r tend to R m * (/) we have 

1 



<7^7T, B<1, 



Rm*{f) Rm*{f) 

which implies that R m *(f) = oo. Let us show that this is not possible. 
In fact, 

[qm* qn,m J qm* Pn,m\ \Z) — -™- n Z "T" " " " , 

and degq m *p n , m < n - 1. It follows that (q m *p n ,m)/qn,m = (q m *Pn,m)/qn,m 
is an incomplete Pade approximant of the function q m *f of type (n,m, 1), 
where the polynomials p n ^ m and q n ,m are relatively prime. As g n ,m(0) ^ 
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0, n > no, the polynomials q m *Pn,m and q n ^ m do not have a common zero at 
z = and \ n = for all n > uq. Again, set m n = degq n , m and 

r n = min {n - 1 - degp njm , m - m n } . 

Notice that r n = m — m n , n > uq, because 

min {n - 1 - deg q m *p n ,m, m - deg g n , m } = 0, n > n . 

Thus, m n +T n = m, n > uq. Using LemmaES] (for m* = 1) and Theorem l2.2l 
we conclude that either Ro(q m * f) < oo or q m * f is a polynomial. However, 
the latter is not possible by hypotheses. On the other hand, Ro(q m *f) < °o 
contradicts R m * (/) = oo. As claimed, / has exactly m* poles in D m * (/). □ 

3. Simultaneous approximation 

Throughout this section, f = (/i, . . . , fd) denotes a system of formal power 
expansions as in ([T]) and m = (mi, . . . , md) € N is a fixed multi-index. We 
are concerned with the simultaneous approximation of f by sequences of 
vector rational functions defined according to Definition 11.11 taking account 
of ([2]). That is, for each n G N, n > |m|, let (R n ,m,i, ■ ■ ■ ,Rn,m,d) be a 
Hermite-Pade approximation of type (n, m) corresponding to f . 

As we mentioned earlier, R nt m t k is an incomplete Pade approximant of 
type (n, |m|, m^) with respect to fk, k = 1, . . . , d. Thus, from (|19p we have 

Dm k Uk) C D\ M (f h ) C D w (/ fc ), k = l,...,d. 

Definition 3.1. A vector f = (fx, ■ ■ ■ , fd) of formal power expansions is 
said to be algebraically independent with respect to m = (mi, . . . , md) £ N 
if there do not exist polynomials pi, ■ ■ ■ ,Pd, at least one of which is non-null, 
such that 

c.l) degp k < m k - 1, k = 1, . . . , d, 

c -2) J2k=iPkfk is a polynomial. 

In particular, algebraic independence implies that for each k = 1, ... ,d, 
fk is not a rational function with at most m^ — 1 poles. Notice that algebraic 
independence may be verified solely in terms of the coefficients of the formal 
Taylor expansions defining the system f . 

Given f = (/i,...,/d) and m = (mi,...,m^) e N d , we consider the 
associated system f of formal power expansions 

T=(fi,...,z m i- 1 f 1 ,f 2 ,...,z m «- 1 fd) = (fi,...J lml ). 

We also define an associated multi-index m given bym=(l,l,...,l) with 
|m| = |m|. The systems f and f share most properties. In particular, poles 
of f and f coincide and R m (f) = R m (f), m e Z + . 

From the definition it readily follows that f is algebraically indepen- 
dent with respect to m if and only if there do not exist constants Ck,k = 
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1, . . . , |m|, not all zero, such that 

|m| 

fc=i 

is a polynomial. That is, f is algebraically independent with respect to m if 
and only if f is algebraically independent with respect to m. By the same 
token, the system poles of f with respect to m are the same as the system 
poles of f with respect to m 

Finally, it is very easy to check that, for all n > |m|, the equations that de- 
fine the common denominator Q n ^ m for (f , m) are the same as those defining 
Qn,m for (f , m) and, consequently, both classes of polynomials coincide. 

Lemma 3.2. Let f = (/i, . . . , fd) be a system of formal Taylor expansions 
as in ([1]) and fix a multi-index m £ N . Suppose that for all n > no the 
polynomial Q n ^m is unique and degQ n ^ m = |m|. Then, the system f is 
algebraically independent with respect to m. 

Proof. Because of what said just before the statement of Lemma [3.2l we can 
assume without loss of generality that m = (l,l,...,l) and d = |m|. We 
argue by contradiction. Suppose that there exist constants Ck, k= l,...,d, 
not all zero, such that X)fc=i c kfk is a polynomial. Should d = 1, Q n ,m = 1 
for all n sufficiently large and deg<5 n ,m < 1 = |m|. If d > 1, without loss of 
generality, we can assume that c\ ^ 0. Then 

d 
fi=P~ ^c fc / fc , 

fe=2 

where p is a polynomial, say of degree N. 

On the other hand, for each n > d — 1, there exist polynomials Q n ,Pn,k, 
k = 2, . . . , d, such that 

- degP n) fc < n- 1, k = 2, ...,d, degQ n <d-l, Q n ^ 0, 

- Qn(z) f k (z) - P n>k (z) = A k z n+1 + • • • , k = 2,...,d. 

Therefore, 

Qn(z)lp(z)-Y J C k f k (z)\ ~ [Qn(z) P (z) -J2 C kPn,k(z)) = Az n+1 + . . . 
\ k =2 J V k=2 J 

and, for n > d + N, the polynomial P n ^ = Q n p — ^2 k= 2 c kPn,k verifies 
degP ni i < n — 1. Thus, for all n sufficiently large, the polynomials P n , k , k = 
1, . . . , d, satisfy Definition 11.11 with respect to f and m. Naturally, Q n gives 
rise to a polynomial Q n ,m with degQ„ im < d = |m| against our assumption 
on Q n ,m- □ 



Set 

Dm(f)=(^| m |(/l)---^|m|(/<i))- 
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The following corollaries are straightforward consequences of Corollary 
12.41 and Theorem 12.61 respectively, together with the fact that, for each 
k = 1, . . . , d, R n ,m,k = Pn,m,k/Qn,m is an incomplete Pade approximant of 
type (n, |m|,mfc) with respect to /&. 

Corollary 3.3. Let f = (/i, . . . , fj) be a system of formal Taylor expansions 
as in ([JJ and fix a multi-index m £ N d . Assume that f is algebraically 
independent with respect to m and there exists a polynomial Q\ m \ of degree 
|m|, Q| m |(0) 7^ 0, such that lim n _ s . 00 Q n<m = Q\ m \. Then Ro(f) > 0, the 
zeros of Q\ m \ contain all the poles that f has in DJ n (f) ; and Ro(fk) < oo for 
each k = 1, . . . ,d. 

Corollary 3.4. Let f = (/i, . . . , f^) be a system of formal Taylor expansions 
as in (H]) and fix a multi-index m = (mi, . . . , m^) € N d . Assume that f is 
algebraically independent with respect to m and there exists a polynomial 
Q\m\ °f degree |m|, Q| m |(0) ^ 0, such that 

limsup ||Q| m | - Q n ,m\\ 1,n = < 1- 

n— >oo 

Then, for each k = l,...,d, either f^ has exactly m^ poles in D mk (fi i .) or 

Ro(Q\m\fk) > Rm k (fk)- 

3.1. Proof of Theorem 11.31 Let us prove first that b) implies a). From 
Lemma 13.21 it follows that f is algebraically independent with respect to m 
and, in turn, from Corollary 13.31 we know that -Ro(f) > 0. So, it is enough 
to prove that f has exactly |m| system poles with respect to m and without 
loss of generality we can assume that m=(l,l,...,l). 

We divide the proof into two parts. First, we collect a set of |m| candidates 
to be system poles of f and prove that they are the zeros of Q\ m \- We also 
prove that any system pole of f must be among these candidates. In the 
second part we prove that all these points previously collected are actually 
system poles of f . 

In the disk Dq({) there cannot be system poles of f since all the functions 
/fe are analytic. Now, for each k = l,...,d, by Corollaries 13.41 and 13. 3[ 
either the disk -Di(/fc) contains exactly one pole of fk, and it is a zero of 
Q| m |, or Ro(Q\ m \fk) > Ri(fk)- Therefore, -Do(f) ^ C and Qi m i contains as 
zeros all the poles of fk on the boundary of -Do(/fc) counting their order for 
k = 1, . . . , d = |m|. Moreover, the functions fk cannot have on the boundary 
of Z?o(/fe) singularities other than poles. 

According to this, the poles of f on the boundary of Dq({) are all zeros of 
Qlml counting multiplicities and the boundary contains no other singularity 
except poles. Let us call them candidate system poles of f and denote them 
by ai, . . . , a ni taking account of their order. Obviously, any system pole of 
f on the boundary of -Do(f) must be one of the candidates since no linear 
combination of the functions in f can produce poles at any other point. 

Since degQi m i = |m| we have n% < |m|. Should n\ = |m| we have found 
all the candidates we were looking for. Let us assume that n\ < |m|. We 
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can find coefficients c±, . . . , C| m | such that 

|m| 

^Ckfk 
k=l 



is analytic in a neighborhood of Dq({). Finding the coefficients Cfc reduces 
to solving a linear homogeneous system of m equations with |m| unknowns. 
In fact, if z = a is a candidate system pole of f with multiplicity r we obtain 
r equations choosing the coefficients c& so that 

(21) I (cj-a) 1 \ f^c k f k (u;) I doj = 0, i = 0,...,r-l. 

j\oj-a\=s y k=1 J 

where 5 is sufficiently small. We do the same with each distinct candidate 
on the boundary of Dq({). The linear homogeneous system of equations 
so obtained has at least |m| — m linearly independent solutions which we 
denote by c 1 -, j = 1, . . . , |m| — n*, n\ < n\. 
Set 

c} = ( c ],i, ■ ■ • > c },|m|)' 3 = !,■■■, I m l - n\. 

Construct the (|m| — nj) x |m| dimensional matrix 



C 1 



r 1 
c l 



C 1 
|m|— n\ 



Define the system gi of |m| — n* functions by means of 

g x = C f = (</!,!,... ,5l,|m|-n*) , 

where (•)* means taking transpose. We have 

|m| 

9i,j = ^2 c },kfk, 3 = 1, ■ ■ ■ , l m l - n\. 
k=i 

As the rows of C 1 are non-null, none of the functions g±j are polynomials 
because of the algebraic independence of f with respect to m = (1, 1, . . . , 1). 
Consider the region 

|m|-raj 

£> (gi)= n ^o^u). 

3=1 

Obviously, by construction, Do(f) is strictly included in £>o(gi) and there 
cannot be system poles of f in Do(gi) \ Do(f). 
It is easy to see that 

|m| 

El r n,m,k 

fc=l ^ n ' m 
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is an (n, |m|,l) incomplete Pade approximant of g\j. Using Theorem 
with m* = 1, for each j = 1, . . . , |m| — n*, either the disk Di(gij) contains 
exactly one pole of gij, and it is a zero of Qi m i, or Ro(Q\m\9i,j) > R±(gij). 
In particular, Do(gi) 7^ C and all the singularities of gi on the boundary 
of I?o(gi) are poles which are zeros of Q\ m \ counting their order. They 
constitute the next layer of candidate system poles of f (now, it is possible 
that some candidates are not poles of f since the functions /& intervene in 
the linear combination as we saw in example ©). All the system poles of 
f on the boundary of L'o(gi) must necessarily be poles of gi for the same 
reason as in the preceding case. 

Let us denote these new candidates by a ni +i, . . . , a ni+ri2 . Of course m + 
n 2 < | m |- Should ri\ + ?i2 = |m|, we are done. Otherwise, ri2 < |m| — n\ < 
|m| — n\ and we can repeat the process. In order to eliminate the n^ poles we 
have |m| — n\ functions which are analytic on Z^o(gi) an d meromorphic on 
a neighborhood of Do(gi)- The corresponding homogeneous linear system 
of equations, similar to (|2ip . has at least |m| — n\ — n<2, linearly independent 
solutions Cj, j = 1, . . . , |m| — n\ — n\, n\ < n%. Set 

—2 ( 2 2 \ • -\ I I * * 

c j = ( c i,i» • • • ' c j,H-n*)> j = l,...,\m\-n 1 -n 2 . 



Construct the (|m| — n\ — n|) X |m| — n\ dimensional matrix 

C 2 = 



r 2 

c l 



r 2 

Define the system g2 of |m| — n\ — ri^ functions by means of 
g2 = Cgl = CCf = (52,1, • • • , 92,\m\-nl-nl) ■ 

The rows of C 2 C X are of the form c 2 C l ,j = 1, . . . , |m| — n\ — n^, where C 1 
has rank |m| — n\ and the vectors c? are linearly independent. Therefore, 
the rows of C 2 C l are linearly independent; in particular, they are non-null. 
Consequently, the components of g2 are not polynomials because of the 
algebraic independence of f with respect tom = (l,l,...,l). Thus, we can 
apply again Theorem 12.61 The proof is completed using finite induction. 

Notice that the numbers m, ri2, ■ ■ ■ which so arise are greater than or equal 
to 1 and on each iteration their sum is less than or equal to |m|. Therefore, in 
a finite number of steps their sum must equal |m|. Consequently, the number 
of candidate system poles of f in some disk, counting their multiplicities, is 
exactly equal to |m| and they are precisely the zeros of Q\ m \ as we wanted 
to prove. 

Now, suppose that there exists a candidate system pole of f that is not 
such or being a system pole has order smaller than the multiplicity of the 
corresponding zero of Q\ m \. Then, for some a S Z+, we have 

Ea = (9a,i,- ■ ■ ,9a,v), v = \m\ -nl - n\ n* , 
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with < n* < rij, j = 0, 1, . . . , a, no = 0, and u = |m|— no— n\ — ■ ■ ■— n a > 
such that there exists a point a on the boundary of the region 

V 

D (ga) = P| D (g aJ ) 

i=i 

that is a pole of order r of g a j for some jo £ {1, ... ,u} but is not a system 
pole of f or is one of order less than r. Let a Ha+ i, . . . , a na+1 , n a+ \ > r, be 
the singularities of the functions g a j on the boundary of Do(g a ) counting 
multiplicities. We distinguish two cases. First, suppose that n\+- ■ -+n a+ i = 
|m|; then n a+ \ = u < v. All the functions g a j admit meromorphic extension 
to a neighborhood of Do(g a )- We pose the problem of finding coefficients 
c\ , . . . , Cj, such that 



E 



Cjdaj 
1 



is analytic on a neighborhood of Do(g a )- The problem consists in solving a 
linear homogeneous system with u equations and v unknowns similar to (|2ip 
but, due to the fact that the point a is not a system pole of f or it is one of 
order less than r, one of the equations may be written as a linear combination 
of the others and we have at most u — 1 equations, with u — 1 < v. So, a 
non-trivial solution necessarily exists which defines a function g analytic on 
a neighborhood of Do(g a ) by means of 

v |m| 

9 = ^2 c i 9a 'i = ^2 dk f k > dk G C ' k = l,...,\m\. 
j=i fc=i 

Following the same argument used in the process carried out to find the 
candidate system poles of f , we deduce that g is not a polynomial. Now, 

|m| 



E 



, -'n,m,fc 

dk 



fc=l 

is an (n, |m|, 1) incomplete Pade approximant of g. Using Theorem 12.61 with 
m* = 1, either the disk D\{g) contains exactly one pole of g, and it is a zero 
°f Qlmh or Ro(Q\m\9) > Ri(d)- But both alternatives are impossible since 
all the zeros of Q\ m \ belong to Do(g a ). So, we have reached a contradiction. 

In case that n\ + • • • + n a +i < |m| we are in the middle of the process 
described above and now, when solving the corresponding system of equa- 
tions to eliminate the n a+ \ poles, we obtain n* +1 < n a+ \ since, again, one 
of the equations is redundant. This implies that, in the last step, say /3, 
when m + • • • + ng = |m| we have |m| — no — n\ — ■ ■ ■ — uq = u < v = 
|m| — n*^ — n\ — ■ ■ ■ — n*n reaching the same contradiction as before. We have 
proved a posteriori that n* = rij, j = 1, 2, . . . 

Thus, the proof of the inverse-type result is complete. Also, we have that 
Q\ m \ = Q|m|(f,m). 
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Let us prove now that a) implies b). Except for some details related to 
the numbers Rt(f,m), where £ is a system pole of f, the arguments are 
similar to those employed in [8j. In spite of this, for completeness, we give 
the entire proof. 

For each n > |m|, let q n<m be the polynomial Q n<m normalized so that 

|m| |m| 

(22) Y^ l-WI = X > Qn, m (z) = ^ K,kZ k - 

fe=l fe=i 

Due to this normalization, the polynomials q ntm are uniformly bounded on 
each compact subset of C. 

Let £ be a system pole of order r of f with respect to m. Consider a 
polynomial combination g\ of type ([6]) that is analytic on a neighborhood of 
Dm except for a simple pole at z = £ and verifies that Ri(gi) = i?^i(f, m) (= 
r^i(f, m)). Then, we have 

|m| 

9i = '^2pk,ifk, degp fc ,i < m k , k = 1, . . . , |m|, 
fe=i 

and 

|m| 

qn,m(z) hi{z) -(*-£) ^2pk,l(z) P n ,m,k( Z ) = ^™ +1 + • ■ ■ > 

k=l 

where h\(z) = (z — £) gi(z). Hence, the function 

q n ,m(z)h 1 (z) z-£ 



z n+1 : 



^lJ2Pk,l(z) P n , m ,ki Z ) 



fc=l 

is analytic on D\{g{). Take < r < R\(gi) and set r r = {z € C : \z\ = r}. 
Using Cauchy's formula, we obtain 



n+l 

du), 



ln,m(u) hi((jj) 



1 f z n+ 
q n , m {z)hi{z)-(z-£) V]pk,i(z)P n!m ,k(z) = — / — — r 

fc=l ■ /lr 

for all z with |z| < r, since deg5^J 1 pfc ) ;|Pn,m,fc < n. In particular, taking 
z = £ in the above formula, we arrive at 

(23) <zn, m (£)M£) = —y^-^——^—^. 

Straightforward calculations lead to 

limSU P |/n(0(Zn,m(0| 1/n <-- 
n— >oo r 

Using that /ii(£) / and making r tend to R\(g\) we obtain 



nmsup|g n , m (0| 1/n < p | , < 1- 
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Now, we employ induction. Suppose that 
(24) limsup ggL(£) 



n— >oo 



Vn |£| 

< o ' frnV J = 0,l,...,s-2 



(recall that i?£ J+ i(f, m) = minjfc=i t ...j+ir^ 1 jfc(f, m)), with s < r, and let us 
prove that formula ()24|) holds for j = s — 1. 

Consider a polynomial combination g s of the type © that is analytic on 
a neighborhood of Dui except for a pole of order s at z = £ and verifies that 
R s (g s ) = *"£,s(f) m )- Then, we have 

|m| 

9s = '^2pk,sfk, degp k)S <m k , fc = l,...,|m|. 
fe=i 

Set h s (z) = (z — £) s g s (z). Reasoning as in the previous case, the function 
q n , m (z)h s (z) z-£ ^ 

z n+l( z _ ty-l z n+l Z^ Pk ' s ^ Z) n ' m ' k{ - Z) 

is analytic on D s (g s ) \ {£}. Put P s = J2k=iPk,sPn,m,k- Fix an arbitrary 
compact set K C {D s (g s ) \ {£,})■ Take 5 > sufficiently small and < 
r < R s (g s ) with K C D r . Using Cauchy's integral formula and the residue 
theorem, for all z € K, we have 

jz)hs(z) 

(z-0 s 



(25) ^rZ-i -t z ~0 Ps(z) = I n {z) - J n (z), 



where 



and 



w-hJ Ir 



z n+l q n , m (u) h s (uj) 

oj n+i (u - zy-i-iw - z) 



T f x 1 /" ^ +1 gn.mMM") , 
J n U = / TT7 Ts T-< r UUJ . 



We have used in (j25|) that deg P s < n. The first integral /„, is estimated as 
in ([23]) to obtain 



(26) limsup \\I n (z)\\r < ^p- = -t^Z-T 

rwoo #«(&) r ^(f, «l) 

As for J n , write 

|m| « (j) rn 

9n,m(w) = 2^ ~\ ( w ~ ■ 

i=o J ' 
Then 

"~ 2 1 /• „n+l „0') 



(27) J n ( z ) = J2^-[ 



z n+l q%U0 M«) ^ 



\ u _ &=9 u n + l j\(u - z) (u - s - 1 - j 
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Using the inductive hypothesis (pit) , estimating the integral in (|2T1) . and 
making e tend to zero, we obtain 



limsup \\Jn{z)\\K < 



\ Z \\K 



\ Z \\K 



lK ~ ^1 ^ >s -i(f,m) ^_i(f,m)' 
which, together with (f26j) and ([25]) . gives 

(28) limsup \\q n ,Uz)hs{z) - {z - if P s (z)\\][ n < J Z } f K 

As the function inside the norm in (|28p is analytic in D s (g s ), inequality 
(J28H also holds for any compact set K C D s {g s ). Besides, we can differentiate 
s — 1 times that function and the inequality still holds true by virtue of 
Cauchy's integral formula. So, taking z = £ in (|28p for the differentiated 
version, we obtain 



lim sup 



\Qn,m'T'!: 



l(«-l) 



(0 



1/n 



< 



ICI 



-R^,s(f,m) 



Using the Leibnitz formula for higher derivatives of a product of two func- 
tions and the induction hypothesis (|24|) . we arrive at 



lim sup 



J--1) 

</n,m 



(0 



1/n 



< 



1*1 



since /i s (£) 7^ 0- This completes the induction. 

Let £i, . . . , £ p be the distinct system poles of f and let r* be the order of 
£i as a system pole, i = 1, . . . ,p. By assumption, n + • • • + r p = |m|. We 
have proved that, for i = 1, . . . ,p and j = 0, 1, . . . , t, — 1, 



(29) 



lim sup 



Qn,m{C,i) 



1/n 



< 



161 



< 



16 



%j+i( f ' m ) ^( f . m ) 



Recall that Qi m i (f , m) is the monic polynomial whose zeros are the system 
poles of f with respect to m. Denote by Lij, i = 1, . . . ,p; j = 0, 1, . . . Tj — 1, 
the fundamental interpolating polynomials at the zeros of Qi m i(f, m); that 
is, for each i = 1, . . . ,p and j = 0, 1, . . . Tj — 1, deg Lij < |m| — 1 and 



■(") 



L\J{b K ) = 5 iK 5 ju , 



l,...,p, v = 0,l,...,n-l. 



Then 

v n-i 

(30) ?n,m(*0 = A n ,| m |Q| m |(f,m) + ^ ^ qn}*n(£i)Lij(z). 

i=\ j=0 

From (gSJ) and flSOJ) it follows that 

lim sup |[g n , m - A n ,| m | Q\ m \ (f , m) ||]/ n < (9 < 1, 

n— >oo 

for any compact if cC, where 



(31) 



= max 



{flJ^r {6P|m|(, ' m) }' 
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As all norms in finite dimensional spaces are equivalent, we obtain 

(32) limsup || gn>m - \n,\m\Q\m\(f,™)\\ 1/n < 9 < 1. 

n—toc 

Now, necessarily we have 

(33) liminf |A n i m i| > 0, 

since if there exists a subsequence A C N such that lim raG A A n i m i = 0, then 
from (J32j) we have lim ne A ||<7n,m|| = 0) contradicting (J22j) . 
As q n>m = \n,\m\Qn,xa, we have proved 

(34) limsup ||Q n , m - Q\ M (f,m)\\V n < 9 < 1, 

n— >oo 

where 6 is given by (f3Tj) . In particular, for n > no, degQn.m = |m|. The 
difference of any two non-collinear solutions Q\ and Q2 of Definition 11.11 
with the same degree and equal leading coefficient produces a new solution 
of smaller degree but we have proved that any solution must have degree 
|m| for all sufficiently large n. Hence, the polynomial Q n ,m is uniquely 
determined for all sufficiently large n. With this we have concluded the 
proof of the direct result. 

Let us prove that the upper bound in (J34|) actually gives the exact rate 
of convergence to obtain ([7]) . To the contrary, suppose that 

(35) limsup ||Q n , m - Q lml (f, m )\\^ n = 9' < 9. 



Let C be a system pole of f such that 

{s|U : « 6 ' w, ' m >}- 



ici „ 1 \(\ 

— — max ■ 



i? c (f,m) {%( 

Naturally, if there is inequality in ([33]) then R((f , m) < 00. 
Choose a polynomial combination 

d 

(36) 9 = ^2pkfk, degp k <m k , k = l,...,d, 

fc=i 

that is analytic on a neighborhood of Dm except for a pole of order s at z = C 
with R s (g) = R( (f , m). On the boundary of D s {g) the function g must have 
a singularity which is not a system pole. In fact, if all the singularities were 
of this type we could find a different polynomial combination g\ of type (|36p 
for which R s {g\) > R s (g) = i?((f,m) against our definition of R^(f, m). 
For short, put Q| m |(f, m) = Q\ m \. Consequently, the function Q\ m \ g can be 
represented as a power series X^o c J zJ w ith radius of convergence R^ (f, m). 
So 



(37) limsup y/\c~\ = l/R c {f, 



m 
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On the other hand, by virtue of ([361) . we have 

d 

H n {z) := Q n , m {z)g{z) - ^Pk{z) P n ,xn,k( z ) = B nz n+1 + ... 

fe=l 
and this function is analytic at least in D\^\ with a zero of multiplicity at 
least n + 1 at z = 0. Taking r < \(\, we obtain 

1 r "»<<"> ^=o. 



2iri J T u n+1 

Set P n = Yjk=\ PkP n , m ,k- Clearly, Q| m | g = (Q| m | - Q n , m ) g + P n + H n and, 
since deg P n < n — 1, we arrive at 

1 /" Q|m|(w)5H 1 /• [Q| m |(w) -Qn.mH] s(w) 

aw = - — - / aoj. 



2m J T uj n+l 2m Jr 

Taking (137p and (|35p into consideration, estimating the integral, and letting 



r tend to |£|, it follows that 



limsup y |c n | < —7 < 



i? c (f,m) „^oo" vl "'-|Cl ICI # C ( f > m ) 
which is absurd. We have completed the proof of Theorem 11.31 □ 

3.2. Convergence of the Hermite-Pade approximants. The following 
result is in some sense the analogue of the formula displayed just after (58) 
in [1] written in different terms. 

Theorem 3.5. Assume that either a) orb) in Theorem \1.3\ takes place. If 
£ is a system pole of order r of f with respect to m, then 



(38) max_ limsup Qn, m (0 

j=0,...,s n ->oo 



i^+i (f , m) ' 
Proof. Let £ be as indicated. From (|29p and (|33p we have 

Vn,m(,c;J 



0,l,...,r-l. 



max_ limsup 

j=0,...,s n->oo 



Vn |£| 

< ff ' fmV 3 = 0,1.. ...r-1. 



J -*-> • 



Assume that there is strict inequality for some s E {0, . . . , r — 1} and fix s. 
Choose a polynomial combination 

d 

g = ^Pkfk, degp k < rn k , k = l,...,d, 

k=i 

that is analytic on a neighborhood of D\%\ except for a pole of order s (< s+1) 
at z = £ with R s (g) = i?£,s+i(f, m). As before, on the boundary of D s (g) 
the function g must have a singularity which is not a system pole. Set 
Qj m |(f, m) = Q| m |. Consequently, the function Qi m i g can be represented as 
a power series X^o c j z ^ with radius of convergence R^g + i(f, m). So 

(39) lim sup vl^J = V^,s+i (f , m) . 
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On the other hand, by virtue of (f36l) . we have 



d 
H n {z) := Qn,m(z)g(z) - ^2pk{z)P n , m ,k{ Z ) = B nZ H+1 + ■■■ 

fc=l 

and this function is analytic in D s (g) \ {£}. Take r smaller than but suf- 
ficiently close to R^a+l (f i m ) an d 5 > sufficiently small. Let Ts, r be the 
positively oriented curve determined by 75 = {co : \co — £| = 5} and T r . We 
have 

1 /■ *■<">.*, = „. 



2vri / ri w™ +1 
Set P n = Yl=iPkPn, m ,k and /i(a;) = (w - £) s sM- Obviously, 

Q|m| 5 = (Q|m| - Qn,m) 9 + Pn + H n 

and, since degP n < n — 1, we obtain 

1 /" QlmlM^M , 1 /" [Qlml ~ Qn,m](u)h{u) 



^n — 


2m 


Jr s , r 


n « w - 

UJ n+1 


1 


JT r 


[Q\m\ - 


Qn,m](w)/i(w) 


2m 


(« 


- £) s W n+1 



27ri7 rir (w-£) s w™ +1 



CtCJ — > - — : / — 4-7 zs n — duj 



^ 2m L, i/!(w - S-I 'a; n + 1 

t/=0 '« 



J_ f [Q\ m \-Qn, m ](u)h(u) y, 1 /■ gSn(6^) 

2vri 7 r? . (a;-^) 8 ^^ 1 W f^2m J ls v\{uj - i) s - u uj n+l 

Estimating these integrals, using (|7|) and the temporary assumption that 



maxlimsup Qn, m (0 

j=0,...,s n _>oo 



v« _. \z\ < \Z\ 



K " i^s+i (f , m) ' 



we obtain 

1 



limsup \c n \ ' n < max < — , 



< 



i%,s+i (f , m) J Rt,s+i (f , m) ' 
which contradicts (f39j) . Hence, (f38j) takes place. □ 



Now, we are ready to give the analogue of (|4|) for simultaneous approx- 
imation. We need to introduce some notation. Fix k € {1, . . . ,d}. Let 
Z)i m j k(f, m) be the largest disk centered at z = in which all the poles of 
fk are system poles of f with respect to m, their order as poles of /& does 
not exceed their order as system poles, and fk has no other singularity. By 
-R|m|,A:(f) m ) we denote the radius of this disk. Let £i, • • • ,£jv be the poles 
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of fk in Di m i fe (f , m). For each j = 1, . . . , N, let Tj be the order of £j as a 
pole of fk and r,- its order as a system pole. By assumption Tj < Tj. Set 

R !m\ fc( f > m ) = min i #|m|,fc(f,m), min i^. ? .(f,m) [• , 

and let D*, fc (f , m) be the disk centered at z = with this radius. 

Recall that o~(B) stands for the 1-dimensional Hausdorff content of the 
set B. We say that a compact set K C C is cr-regular if for each zq € K and 
for each 5 > 0, it holds that a{z € K : \z — zq\ < 5} > 0. 

Theorem 3.6. Let f be a system of formal Taylor expansions as in ([1]) and 
fix a multi-index m € Z+ \ {0}. Suppose that either a) or b) in Theorem 
\1.3\ takes place. Then, 

(40) limsup ||/ fc - -R w ,m,fc||^ n < „« Tyr , fe=l,...,d, 

n->oo "*m k \Jk) 

where K is any compact subset of DTAf^) \"P| m i(f, m). //, additionally, K 
is a -regular, then we have equality in (j40|) . Moreover, 

< ]b (/fc) = -Rfm|,fc(f,m)» fc = l,.--,d- 

Proof. Let us fix A; € {1, . . . , a!} and maintain the notation introduced above. 
Let K be a compact subset contained in -Dj* m i fe (f, m) \ "P| m |(f, m). Take r 
smaller than but sufficiently close to iff. fe (f, m), and (5 > sufficiently 
small so that K is in the region bounded by T r and the circles {z : |z — £j| = 
(5}, jf = 1, . . . , N. Let Ts t r be the curve with positive orientation determined 
by T r and those circles. On account ofDefinition ll.il using Cauchy's integral 
formula we have 

{Qn,mJk ~ Pn,m,k){z) = T— / — — T « w 

Since lira n Qn,m = Q\ m \i using (j38j) and standard arguments we obtain 

(41) limsup || / fc - Rn,ia,k II x ^ 



«f m l, fc ( f > m )' 

This last relation implies that 0"-lim n _ s . oo -R n>m> fc = fk inside DT, fc (f , m). 
Since Rm k (fk) is the largest disk inside of which such convergence takes place 
it readily follows that RT, fc (f, m) < R* mk (f k ). Should DT, fc (f,m) contain 
on its boundary some singularity which is not a system pole then necessarily 
-Rfi fc (f, m) = Rm k {fk) because a-convergence implies that all singularities 
inside must be zeros of Q\ m \ but the zeros of this polynomial are all system 
poles as we proved in Theorem 11.31 Assume that Rm k (fk) > R^mi fc(f> m )- 
Then, we have R^Afk) > m h 1 i=l ) ...,iV-%yr(f> m )- From the proof of [2J 
Theorem 3.6] we know that for each pole £ of order r of /& inside D*{f) 



k 



limsup Qn, m (0 



V« |6| 
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This contradicts (|38|) . Consequently Rm k (fk) = -Rf m i fc(f> m ) as claimed. 

Due to (|4ip , we have also proved (|40p . In order to show that this formula is 
exact for <r-regular compact subsets one must argue as in the corresponding 
part of the proof of [2j Theorem 4.4]. □ 



As compared with [2J Theorem 4.4], Theorem 13.61 offers weaker assump- 
tions and a characterization of the values R*{fk) in terms of the analytic 
properties of the functions in the system instead of the coefficients of their 
Taylor expansion. An open question is to obtain an analogous characteriza- 
tion when the assumptions of Theorem 13.61 do not take place. 

It would be interesting to study inverse problems for row sequences of 
Hermite-Pade approximation when only the limit behavior of some of the 
zeros of the polynomials Q n ,m is known, in the spirit of the conjectures 
proposed by A. A. Gonchar in pQ. 
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